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Summary 

^"^The  pressure  and  displacements,  due  to  line-force  and  line-source 
excitation,  aze  given  in  terms  of  Pourier  integrals  Which  can  be  evaluated 
numerically.  Plots  of  intensity  vectors  provide  a  vivid  illustration 
of  the  energy  flow  in  the  fluid.  The  selected  plots  show  significant 
interchange  of  propagating  energy  between  fluid  and  the  driven  plate  when 
the  excitation  is  a  force,  but  not  when  the  excitation  is  a  source. 
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In  the  absence  of  dissipation,  the  total  energy  flow  along  a  uniform 
pipe  containing  a  stationary  fluid  is  conserved.  However,  because  the 
system  normal  modes  are  non— orthogonal  over  the  fluid's  cross-section,  the 
distribution  of  energy  between  the  pipe's  wall  and  the  fluid  need  not  be  a 
constant  with  distance  along  the  pipe.  This  interchange  of  energy  between  a 
pipe  and  its  fluid  is  not  without  interest  to  the  noise  control  engineer 
who  requires  to  understand  the  physics  of  energy  interchange  between  mech¬ 
anical  and  fluid  elements  of  a  dynamic  system.  James  [1]  has  given  a  qual¬ 
itative  description  of  the  energy  flow  in  a  uniform  pipe,  with  time— harmonic 
point  source  and  force  excitation,  in  terms  of  intensity  vector  plots  in  the 
fluid  which  show  clearly  this  energy  interchange.  Puller  [2]  has  provided  a 
quantitative  description  of  this  phenomenon  by  giving  the  mathematics,  and 
some  numerical  results,  of  the  energy  distribution  between  pipe  and  fluid. 

The  purpose  of  the  work  described  herein  is  to  give  the  preliminary 
work  necessary  for  a  combined  analytical  and  numerical  investigation  of  the 
simpler  problem  of  the  time-harmonic  line-force  or  line-source  excitation  of 
a  system  comprising  of  a  layer  of  fluid  bounded  by  two  identical  thin  plates 
This  problem  is  much  more  tractable  than  the  corresponding  pipe  problem,  and 
its  greater  simplicity  may  help  to  clarify  the  physical  processes  that  are 
involved  in  the  energy  interchange  between  fluid  and  plates.  This  initial 
work  will  give  the  basic  mathematics  needed,  and  also  discuss  some  plots  of 
acoustic  intensity  vectors  which  provide  a  vivid  illustration  of  the  energy 
flow  in  the  fluid. 


2.  MATHEMATICS  OP  PROBLEM 

The  geometry  of  the  problem  is  shown  in  Figure  1.  A  time-harmonic 
line-force  or  line-source  drives  a  system  comprising  of  a  layer  of  fluid 
contained  between  two  identical  infinite  elastic  plates.  The  differential 
equations  satisfied  by  the  plate  displacements  and  fluid  pressure  are 


[D a4/ax4  -  w2p  h]*r<x)  «  F  8(x)  +  p(x,H)  (1) 

8  1  1 

[D34/9X4  -  «2psh]M2(x)  -  F28(X)  -  P(X,0)  (2). 

{92/9z2  +  92/ dx2  +  k2]p<x,Z)  -  -4trp08(x>6(z-z0)  (3) 


where  D-Eh3/12( 1-c2 )  is  the  flexural  rigidity  of  the  plates  whose  densities 
are  p8  and  whose  thicknesses  are  h;  W^x)  and  tf2(x)  are»  respectively,  the 
vertical  displacements  of  the  upper  and  lower  plates;  p(x,z)  is  the  acoustic 
pressure  in  the  fluid  layer  whose  density  is  p  and  whose  sound  velocity  is 
c;  k  is  the  acoustic  wavenumber,  u/c;  Ft  and  f2  are  the  amplitudes  of  the 
line-forces  that  are  applied  to  the  upper  and  lower  plates,  respectively; 
p0  is  the  amplitude  of  the  line-source  whose  free-field  pressure  is  given 
in  numerous  texts,  for  example  [3],  as 


Gp(x,z)  •  ffip0H0{kv(x  +z-zQ  )> 


(4) 


In  which  H0  is  the  Hankel  function  J0+iY0>  the  time-harmonic  factor 
exp( -iut ) ,  is  omitted  throughout. 


The  analysis  proceeds  via  the  Fourier  transform 


where 
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(6) 


The  general  solution  of  the  differential  equation  (3)  in  the  spectral 
domain  is 


where  [3] 


and 


p(a,z)  -  Gp(«,z)  +  P8(a.z) 

Gp(a,z)  -  2»riP0y-lexp(iy(z-z0)>,  z>zQ 
Gp(a,s)  -  2trip0y”1ejcp{-iy(*-*0)},  z<zQ 

Ps(a.«)  -  A<a)exp(iyz)  ♦  B< a >exp( -iyz ) 


(7) 


(8) 


in  which  y-+v(k2-a2),  and  A(a)  and  B(a)  are  the  unknown  constants  of  integ¬ 
ration  of  the  homogeneous  wave-equation. 

the  spectral  forms  of  equations  (1)  and  (2)  may  be  written  in  matrix 
notation  as 


(Do’  -  u‘pgh)  W1(o)  -  Fj  +  P<a,H) 


W2(a) 


-p(a,0) 


Eliminate  p  from  this  equation ,  by  the  use  of  equations  ( 7 )  and  ( 8 )  to  give 


(Da4  -  «2pgh)  l?1(a)  -  F1  +  exp(iaH)  exp(-ioH)  A(a) 

W2(a)  P2  -1  -1  B(a> 


+  2  rrip^y  exp  {  iy(  tt-zQ  ) } 


-exp(iyzQ) 


The  boundary  conditions  at  the  interfaces  require  continuity  of  normal 
displacements,  viz. 


{dp (x,z)/dz> 


{dp(x,z)/dz}z_0  -  pofa^x) 


which  give  the  spectral  equations 


»i(o)  -  iy  exp( iyH)  -exp(-iyH)  A(a)  +  2trp0  -exp{iy(tt-z0)} 


W2(<*) 


exp( iyzQ ) 


Elimination  of  wt( a)  and  W2( a )  between  equations  (10)  and  (12)  gives 
the  matrix  equation  from  which  the  unknowns,  A (a)  and  B(a),  may  be  obtained, 
viz. , 


iy  D+exp(iyH)  -D_exp(-iyH)  A(a)  -  pw2  +  2trpQD+  exp{iy(H-zQ)) 
D~  -D+  B(a)  P2  -exp(  iyzQ ) 


where 


D  -  Da  -  (i)  p  h  +  ipw  /y 


D  «  Da4  -  «2p#h  -  ipw2/y 


n»  spectral  pressure,  p(a,z),  is  obtained  from  equation  (7),  and  the  spectral 
displaceemnts,  w1(a)  and  W2( a)  are  obtained  frost  equation  (12).  The  spectral 
acoustic  particle  displacements,  Wx( a, z )  and  wz(a,z),  in  the  horizontal  and 
vertical  directions,  respectively,  are  given  by  the  equations 


where 


t*>\(  a.  z ) 


p»2Wz(a,z) 


C(z,zQ) 


iap(  a,  z ) 


iy(A( a)exp( iyz )  -  B< a )exp<-iyz ))  +  C(z,zQ) 


-2trp0eiq>(iy(z-z0)},  z>zQ 


C(z,zQ)  »  2tip0e*p{-iy(z-z0)},  z<zQ 


The  displacements  of  the  plates,  and  the  pressure  and  particle  displacements 
in  the  fluid,  are  obtained  by  numerical  or  contour  Integration  of  the  Fourier 
integrals . 


The  Fourier  integrals  of  equation  (4),  with  upper-  limits  of  integ¬ 
ration  set  to  approximately  twice  the  highest  free-wavenumber  at  the 
selected  frequency,  were  evaluated  by  a  simple  adaptive  Gaussian  quadrature 
scheme.  Because  the  integrals  must  be  evaluated  in  the  principal  value 
sense,  it  is  necessary  to  introduce  damping  into  the  system  via  a  complex 
Young's  modulus,  E(  l-iT)g )»  and  a  complex  sound  velocity,  c(l-ihf).  The 
following  constants  in  SI  units  were  used  in  the  computations: 

E«19.5xl010  ff-O. 29  ps-7700.0  h-0.01  or  0.05 

p-1000.0  C— 1500 . 0  H— 0.20 
r)8«0 . 02  T)f"0  •  001 

The  Chosen  values  of  the  loss-factors  have  little  effect  on  the  appearance 
of  the  intensity  vector  plots. 


Some  dispersion  plots  of  the  dissipation- free  system  are  shown  in 
Figures  2  and  3  for  the  cases  of  plate  thickness  of  0.01  and  0.05m,  resp¬ 
ectively.  In  Figure  2,  the  near  identical  branches,  labelled  1  and  2,  are 
symmetric  and  antisymmetric  modes  Whose  motion  is  mostly  confined  to  the 
plates,  their  group  velocity  being  approximately  twice  their  phase  velocity. 
The  branches  labelled  3-5  are  predominantly  fluid  waves  in  which  the  phase 


and  group  velocities  at  cut-on  are  infinity  and  zero,  respectively:  in  fact, 
the  branch  labelled  3,  at  frequencies  above  3kHz,  is  close  to  a  plane-wave. 

In  Figure  3,  the  branches  labelled  1  and  2  start  as  predominant ly  plate 
leaves  in  which  the  group  velocity  is  approximately  twice  the  phase  velocity, 
and  end  as  fluid  waves  in  which  the  phase  and  group  velocities  are  approx¬ 
imately  equal.  The  branches  labelled  3-5  start  as  fluid  waves,  but  they  also 
Change  their  nature  as  the  frequency  increases.  It  Should  be  kept  in  mind, 
however,  that  the  plate  theory,  for  the  0.05m  plate,  is  of  limited  validity  at 
frequencies  above  5kHz.  Crighton  [4]  has  discussed  the  free-waves  on  a  fluid 
loaded  elastic  plate  and  Fuller  &  Fahy  [5]  have  discussed  the  changing  nature 
of  the  wavenumber  branches  with  frequency  for  the  case  of  a  fluid-filled  pipe. 


The  components  of  the  intensity  vectors  in  the  fluid-layer  sure  defined 


as 


•  * 


Ix(x,z)  -  (1/2)  Real(p(x,z)Wx(x,z)] 


(17) 


I  (x,z)  -  (1/2)  Keal[p(x,z)W  (x,z)] 
z  z 


where  *  denotes  complex  conjugate.  The  vector  plots  at  selected  frequencies 
are  shown  in  Figures  4-11  in  which  the  even  numbered  Figures  are  for  force 
excitation  of  the  lower  plate,  while  the  odd  numbered  Figures  are  for  a 
source  excitation  located  at  r-H/3 .  Due  to  the  symmetry  of  the  problem, 
only  the  vectors  for  x>0  are  shown;  they  are  normalised  to  the  same  maximum 
length  in  each  Figure,  their  individual  lengths  being  proportional  to  I1/2. 
Vectors  are  plotted  only  for  a  plate  thickness  of  o.oim;  hence,  it  is  Fig.  2 
which  shows  the  free -wavenumbers  of  the  system. 

In  Figure  4,  the  interaction  of  the  symmetric  and  antisymmetric  plate 
waves  (1,2)  produces  propagating  energy  flow,  at  distances  in  excess  of  0.4m, 
which  decays  exponentially  in  the  vertical  direction.  The  circulating  energy 
flow  at  0.2m  is  caused  by  interaction  of  propagating  modes  with  evanescent 
modes;  they  are  a  common  characteristic  in  intensity  vector  plots.  In 
Figure  5,  at  distances  less  than  0.4m,  the  absence  of  propagating  fluid  waves 
results  in  a  dominance  of  non-propagating  energy  flow  caused  by  interaction 
among  evanescent  modes;  the  weak  excitation  of  the  plate-waves  (1,2)  being 
evident  at  distances  in  excess  of  0.4m. 

In  Figure  6,  the  interaction  between  the  plate-waves  (1,2)  and  the 
fluid-wave  ( 3 )  results  in  a  well-defined  interchange  of  energy  between  the 
fluid  and  the  lower  plate:  also  evident  in  the  upper-half  of  the  layer  is 
the  fluid-wave  (3),  on  its  own,  which  is  close  to  a  plane-wave.  In  Figure  7 
the  fluid-wave  ( 3 )  dominates  the  energy  flow,  the  interaction  between  plate 
and  fluid  being  insignificant. 

The  increasing  complexity,  with  increasing  frequency,  of  the  energy 
flow  is  shown  in  Figures  8-11.  Of  particular  interest  is  the  rising  and 
falling  of  the  energy  in  the  fluid;  the  well-defined  interchange  of  energy 
between  fluid  and  plate  when  the  excitation  is  a  line-force;  and,  the 
absence  of  significant  energy  interchange  when  the  excitation  is  a  source. 


Sons  vector  plots  have  also  been  obtained  for  the  case  of  the  fluid 
layer  bounded  by  0.05a  plates,  but  they  are  not  shown  here.  At  10kHz, 

Which  is  above  the  plate  ‘coincidence*  frequency,  there  is  a  very  strong 
interchange  of  energy  between  the  upper  and  lower  plates.  These,  and  other, 
intensity  vector  plots  will  be  shown  in  continuation  reports. 


4.  COMCUJPIHG  REMARKS 

Formulae  have  been  given  for  the  Fourier  transforms  of  the  pressure 
and  displacements,  due  to  tiae-harmonic  line-force  and  line-source  excitation. 
Some  plots  of  intensity  vectors  show  propagating  energy  interchange  between 
the  fluid  and  the  force  excited  plate,  which  is  a  result  of  an  interaction 
between  propagating  plate  and  fluid  waves.  Vectors  due  to  source  excitation 
do  not  show  energy  interchange,  due,  presumably,  to  the  very  weak  excitation 
of  the  plate  waves.  The  plots  provide  a  vivid  illustration  of  the  energy 
flow  in  the  fluid,  their  physical  interpretation  being  greatly  facilitated  by 
the  availability  of  wavenumber— frequency  plots. 

FOllow-up  projects  should  include;  (a)  further,  but  trivial,  analytical 
work  to  try  to  obtain  simple  closed-fora  expressions  for  the  transformed 
quantities  of  interest,  for  example,  it  is  possible  to  obtain  the  dispersion 
relation  from  the  determinant  of  equation  (13)  as 


tan(>H)/y  +  (Da4  -  w2pah)/{y2(0a4  -  «2p> h)2  -  p2u>4}  -  0;  (18) 

(b)  numerical  evaluation  of  the  energy  distribution  between  fluid  and  plates, 
as  a  function  of  horizontal  distance,  for  a  variety  of  material  and  geometric 
constants;  (c)  interpretation  of  mathematics  and  numerical  results  to  give 
conditions  under  Which  energy  transfer  between  plates  and  fluid  is  of  most 
significance;  (d)  investigation  of  the  effects  of  line-constraints ,  which 
can  be  attached  to  the  plates  by  the  method  of  dynamic  stiffness  coupling 
[6] i  and  finally,  (e)  extension  of  the  work  to  include  the  effects  of 
plates  of  unequal  thicknesses,  with  fluids  rather  than  vacuums  in  the 
upper  and  lower  half-spaces,  such  a  project,  based  partly  on  Fuller's 
paper  [2]  and  the  work  herein,  may  help  considerably  to  clarify  the  physics 
of  energy  flow  in  fluid-filled  pipes. 


J.H.  James  (FSO) 
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